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Starting from a spin-orbital model for doped manganites, we investigate a competition between ferromag- 
netic and antiferromagnetic order in a one-dimensional model at finite temperature. The magnetic and 
orbital order at half filling support each other and depend on a small antiferromagnetic superexchange be- 
tween t2 g spins and on an alternating Jahn-Teller potential. The crossover to a metallic ferromagnetic phase 
found at finite doping is partly suppressed by the Jahn-Teller potential which may localize e g electrons, 
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1 Introduction The transition from an insulating to metallic behavior in doped manganites, observed 
as a function of temperature or filling, is one of the outstanding problems in the field of strongly correlated 
systems. The theoretical challenge is to understand a rich variety of metallic, insulating, magnetically 
ordered and orbital ordered phases of doped manganites in terms of the dynamics of correlated e g electrons 
which involves their orbital degrees of freedom 1 1 1. Magnetic interactions: spin superexchange (SE) and 
double exchange (DE), as well as orbital SE, frustrate each other and their competition leads to several 
magnetic phases, including phase separation I2l . 

While orbital order supports A-type antiferromagnetic (AF) phase in LaMnC>3 0, one expects that 
e g orbitals play an important role in the insulating ferromagnetic (FM) phase at finite doping, below the 
crossover to a metallic FM phase 1 1 1. Here we investigate a one-dimensional (ID) model derived from the 
realistic spin-orbital model for manganites 1 3 1, showing an interplay between spin and orbital order. 

2 Effective orbital t-J model at finite temperature We investigate the ID orbital t-J model for 
manganites for a given configuration of core spins S, 

H{S)=H t + Hj + Hj,+H n , (1) 

which consists of the hopping term H t , the orbital superexchange Hj for the e g electrons, superexchange 
Hji for the core spins formed by t2 g electrons, and an alternating Jahn-Teller (JT) potential Hjt- 

As we consider ID chains in z-direction, only electrons within \z) = -j=\3z 2 — r 2 ) orbitals are mobile, 
while those within \x) = -^\x 2 — y 2 ) orbitals are localized: 

Ht = -^2u,i+i{4z 5 i+i,z + 4+i.Az), (2) 

i 

where Ci Z and c] are the annihilation and creation operators restricted to the Hilbert space without double 
occupancies. Due to strong Hund's coupling between e g and t2 g electrons, the hopping amplitude, Uj+i = 
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tui^+i = t cos(0i ! i+i/2)e* Xii+1 , depends of the relative angle 6i ^+1 of the two classical core spins at sites 
i and i + see e.g. Ref. |4|. The complex phase Xi.i+i can be neglected for a ID chain |5 1. 

It is convenient to introduce orbital pseudospin operators Tf = \a\ = \{ni X — rii Z ), with eigen- 
states \z) and \x). With the electron number operator fii restricted to the Hilbert space without double 
occupancies, the superexchange term which follows from Refs. 1 3 1 and 1 6 1 becomes 

Hj = ( 2 \ u \li+i + Z)( 2T i T i+i - \^n l+1 ) - ^JJ2 C 1 - Hli+i)nizn i+ i, z 



i 1 - Mm+i) n iz {l-n i+ i) + (l-ni)n i+ i, z j. 



(3) 



Like the kinetic energy, it depends on the direction of core spins via the factors \u\f The first term 
oc (2\u\f i+1 + 3) favors FM spin order for ion pairs Mn 3+ -Mn 3+ , while the second and third ones oc 
(1 — |u|f favor AF order for Mn 3+ -Mn 3+ pairs [oc n iz n i+ i_ z ], and for Mn 3+ -Mn 4+ pairs [oc n iz {l— 
ni+i)]- The SE parameter J = t 2 /e( 6 Ai) is given by the high-spin excitation energy e( 6 Ai) and the 
prefactors of the above terms were chosen to preserve the relative magnitude of the different excitations for 
realistic parameter values [3 1. Throughout this work, we use J = 0.1252; and choose t as unit of energy. 
The third term in Eq. (1) is an AF Heisenberg interaction between the core ti g spins, 

Hj, =J'Y / (S l -S t+1 -S 2 ). (4) 

i 

As we used classical core spins Si of unit length, their physical value 5 = 3/2 was compensated by a 
proper increase of J'. Finally, we include an alternating JT potential of the form 

H n = 2E n eM^Ri) T*, (5) 

i 

We employ a Markov chain Monte Carlo (MCMC) algorithm for the classical core spins (see e.g. Ref. 
|4|) combined with exact diagonalization for the electron degrees of freedom |7|. In this method, the 
many -particle Hamiltonian determined by {u^.i+i}, which are in turn given by the core spin configuration 
S, is solved by exact diagonalization taking into account its block-diagonal structure. The lowest states of 
each block are used to evaluate the trace over the fermionic degrees of freedom, Tr c e~' 3W (' s ) —: w(S), 
which gives the statistical weight for S and which is sampled by the MCMC. Autocorrelation analysis was 
used in order to ensure that enough configurations were skipped between measurements. 



3 Numerical results and discussion We begin with analyzing spin and orbital correlations at increas- 
ing doping. Figure^shows the core spin structure factor, S(k) = Y^ij SiSj£~ lk ( l ~i\ for an N = 12 
chain and doping 0<x<5/12(j: = l — n, where n is electron filling) near the system filled by n = 1 
electron per site (as in LaMnC»3), and for three different parameter sets at low temperature. The orbital 
structure factor, T(k) = ■ T/TJe - ^^ - - 7 ^ /n 2 , which is renormalized by the squared electron density 
for better comparability of results obtained at different fillings, for the same parameters is shown in Fig. 
[2] Altogether, the orbital and magnetic order obtained at x = and Ejt = complement each other — 
polarized (ferro) \z) orbitals at J' = 0.02t induce AF spin order [Figs. [Qa) and|2ja)], but if the occupied 
e g orbitals alternate, FM spin order is supported instead [Figs, ^b) and|3b)]. Increasing doping induces 
a gradual crossover towards FM order due to the DE mechanism, being the most efficient in the metallic 
phase when only \z) orbitals are occupied. This insulator-metal transition is hindered by the JT potential 
which favors orbital alternation [Fig. |2jc)] and suppresses FM correlations [Fig. [3c)] . 

For a finite SE between core spins J' = 0.02i and in absence of JT potential (Ejt = 0), a transition 
from AF to FM spin order is induced by doping [Fig. |l(a)| , and predominantly the mobile \z) orbitals 
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Fig. 1 Core spin structure factor S(k) for: (a) J' = 0.02t, E JT = 0; (b) J' = 0, £j T = 0; (c) J' = 0.02t, E JT = 
O.lt, and several doping levels, as obtained for an N = 12 site chain. Parameters: J = 0.125t and fit = 100. 
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Fig. 2 Orbital structure factor T(k) for: (a) J' = 0.02t, Ejt = 0; (b) J' = 0, £.j T = 0; (c) J' = 0.02t, 
Ejt = O.lt, and several doping levels, as obtained for an TV = 12 site chain. Parameters: J — 0.125t and fit — 100. 



are occupied at all doping levels [see Fig. |2(a)j . In the spectral density depicted in Fig. |3(a)| one sees 
signals with low weight for the \x) orbitals and a rather broad incoherent peak due to the \z) excitations. 
The broad peak contains dispersionless states stemming from \z) electrons next to the few |x) electrons. 
Additionally, a dispersive band is found which comes from \z) electrons moving in the AF background, 
showing a remarkable similarity to a one-orbital model [5 1. At the higher temperature j3t — 30, more |sc) 
orbitals are occupied (4.53 ± 0.07 ~ 38% vs. 2.03 ± 0.1 at fit = 100), and their weight in the spectral 
density increases [Fig. |3(b)| . The dispersionless broad feature becomes stronger, as the motion of \z) 
electrons is hindered by \x) electrons. Both effects are further enhanced when spin order weakens at very 
high temperature f3t = 10. 

A chain doped with three holes x = 0.25 has FM order of core spins [see Fig. |l(a)| and practically 
only mobile \z) electrons [Fig. 12(a)] , due to the dominating DE mechanism which favors both FM order 
and ferro orbital polarization. The spectral density at low temperature jit = 100 is shown in Fig. |3(d)| and 
one sees a coherent band with almost unrenormalized dispersion ~ At due to e g electrons moving in \z) 
orbitals, and only a weak signal above the Fermi energy for the immobile excitations in \x) orbitals. For 
j3t — 30 [see Fig. |3(e)| , the peaks of the |,z}-band are broadened, the bandwidth is slightly reduced and 
a small spectral weight is induced in \x) states below the Fermi energy (0.12 ± 0.02 for the filling by 9 
electrons). Again, these effects become more pronounced for the high temperatures j3t — 10, see Fig. |3(f)| 
where 1.18 ± 0.06 electrons occupy \x) orbitals. The broadening of the peaks and the reduction of the 
bandwidth result from core spin fluctuations weakening FM correlations, while the additional (incoherent) 
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Fig. 3 Spectral functions A(k,uj) for an N = 12 orbital chain (1). First row — undoped chain (x = 0) at tem- 
peratures: (a) fit = 100, (b) pt = 30, and (c) f3t = 10; second row — doped chain with three holes (x — 0.25) at 
temperatures: (d) fit = 100, (e) /3t = 30, and (f) /3t — 10. Solid and dashed lines for \z) and \x) excitations. Insets 
show the density of states n(u). Parameters: J = 0.125t, Ejt = 0, J = 0.02i. 



signals besides the tight-binding band are caused by the \x) electrons. The reduction of FM correlations 
and the increasing electron density within |x) orbitals occur simultaneously with rising temperature. 

Yet another situation arises (at J' = 0.024) in presence of a small alternating JT potential Ejt = O.li. 
It favors alternating orbital pattern [Fig. |2(c)| , which in turn destroys the core spin AF order in the undoped 
chain [Fig. 1 1(c)] , Upon doping, the initially weak FM correlations increase, \x) orbitals are depleted and 
the number of\z) electrons even increases slightly to ~ 7 for 2 and 3 holes. For the undoped chain there are 
just three peaks in the spectral function A(k, u>) [Fig. |4(a)| : a central peak from the \z) electrons trapped 
between two occupied \x) orbitals, and two side peaks due to the \x) electrons. Very similar results (not 
shown) are obtained for J' = and Ejj = 0, where the orbitals also alternate, see above. 

Doping with one hole removes an \x) electron and a three-site well is created, leading to three \z) signals 
in the spectral density depicted in Fig. |4(b)| The spectral function A(k, uj) for doping with three holes [Fig. 
|4(c)| shows new features on top of the tight-binding-like band at Ejt = (Fig. |3(d)| >, that arise from the 
two remaining \x) electrons and from the alternating on-site potential seen by \z) electrons. Note that the 
FM is weaker at this doping level than when the JT potential is absent, because DE is hindered. From 
doping x > 1/3 onward, the chain is fully polarized and FM order is enhanced, see Figs. |l(c)| and |2(c)| A 
strong JT potential Ejj ~ 0.5< delays a crossover to the metallic FM phase to very high doping x > 0.4. 
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Fig. 4 Spectral functions A{k,uS) for an iV = 12 orbital chain as in Fig. [5] but with JT potential En — O.li, as 
obtained for: (a) undoped chain (x — 0), and chain doped with (b) one hole x ~ 0.083, and (c) three holes x = 0.25. 

Finally, we determined the DE constant for the present ID chain, Jde = (H t ) / (2zS 2 ), where z = 2 
and S = (4 — x)/2 |6|. While the DE increases with x and is simply frustrated by the AF SE between 
core spins (Jde does not depend on J'), increasing JT potential reduces Jde in the insulating regime. For 
instance, at x = 0.25 one finds J DE ^ 0.030, 0.027 and 0.02U for £? JT = 0, 0.1 and 0.25i However, 
at doping x = 0.42 the FM phase is already metallic for the entire range of < Ejt < 0.25i, with 
J DE ^ 0.046*. 

4 Summary and conclusions Even in one dimension, the orbital t-J model shows an interesting com- 
petition of FM and AF SE terms. At low temperatures and for the undoped chain, a relatively small AF SE 
J' between the <2g core spins already suffices for a transition from FM spin order and alternating orbital 
order (J' = 0) to AF spin order and ferro orbital polarization (J' = 0.02i). These different situations, both 
encountered qualitatively in LaMnC>3, result in quite distinct one-particle spectra. Because of the close 
interplay of orbital- and spin-degrees of freedom, the change between these two scenarios can also be trig- 
gered by an alternating JT potential (Ejt — O.li). At higher temperatures, the competing interactions lead 
to further changes in the spectra as both the spin and the orbital order are destroyed. 

For the doped chain (x = 0.25) without JT potential, the kinetic energy favors FM order and orbital 
polarization. This metallic FM phase may be destroyed either by finite JT potential, or by increasing 
temperature, when both the FM order and the orbital polarization are reduced, leading in both cases to in- 
coherent spectra with reduced bandwidth. This behavior indicates that an insulating phase may be induced 
either by increasing temperature, or by enhancing the JT distortions by chemical substitution. Summa- 
rizing, the present study highlights: (i) a complex interplay between spin, orbital and change degrees of 
freedom, and (ii) the importance of the coupling to the lattice — all these all these aspects decide about 
spin and orbital correlations in doped manganites and have to be included in their realistic microscopic 
models|8|. 
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